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ABSTRACT 

We explore the relation between the structure and mass accretion histories of dark 
matter halos using a suite of cosmological simulations. We confirm that the formation 
time, defined as the time when the virial mass of the main progenitor equals the mass 
enclosed within the scale radius, correlates strongly with concentration. We provide a 
semi-analytic model for halo mass history that combines analytic relations with fits 
to simulations. This model has the functional form, M{z) = Mq{1 + where 

the parameters a and j3 are directly correlated with concentration. We then combine 
this model for the halo mass history with the analytic relations between a, /3 and the 
linear power spectrum derived by Correa et al. to establish the physical link between 
halo concentration and the initial density perturbation field. Finally, we provide fit¬ 
ting formulae for the halo mass history as well as numerical routineflj, we derive the 
accretion rate as a function of halo mass, and demonstrate how the halo mass history 
depends on cosmology and the adopted definition of halo mass. 
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1 INTRODUCTION 

Dark matter halos provide the potential wells inside which 
galaxies form. As a result, understanding their basic proper¬ 
ties, including their formation history and internal structure, 
is an important step for understanding galaxy evolution. It 
is generally believed that the halo mass accretion history 
determines dark matt er halo properties, s uch as their ‘uni¬ 
versal’ density profile (iNavarro et al.|[l996l . hereafter NFW). 
The argument is as follows. During hierarchical growth, ha¬ 
los form through mergers with smaller structures and accre¬ 
tion from the intergalactic medium. Most mergers are minor 
(with smaller satellite halos) and do not alter the structure 
of the inner halo. Major mergers (mergers between halos 
of comparable mass) can bring material to the centre, but 
they are found not to play a pivotal role in mo difying the 
internal mass distribution (IWang fc Whitell2009l l. Halo for¬ 
mation can therefore be described as an ‘inside out’ process, 
where a strongly bound core collapses, followed by the grad¬ 
ual addition of material at the cosmological accretion rate. 
Through this process, halos acquire a nearly universal den- 
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sity profile that can be described by a simple formula known 
as the ‘NFW profile’ (NFW). 

The origin of this universal density profile is not fully 
understood. One possibility is that the NFW profile re¬ 
sults from a relaxation mechanism that produces equilib¬ 
rium and is largely independent of th e initial conditions 
and merger history dNavarro et al.l[l99itl . Ho wever, another 
popul ar explanation, originally proposed by ISver fc White! 
(Il998ll . is that the NFW profile is determined by the halo 
mass history, and it is then expected tha t halos should also 


follow a universal mass history profile ( 

Dekel et al.l 20031: 

Manriaue et al.l l2003l: ISheth & TormenI 

2004 

Dalai et al.l 

2010l: Salvador-Sole et al. 20121: Giocoli et al. 

I2OI2II. This 


universal accretion history was recently illustrated by 
iLudlow et al.l (l2013f l. who showed that the halo mass histo¬ 
ries, if scaled to certain values, follow the NFW profile. This 
was done by comparing the mass accretion history, expressed 
in terms of the critical density of the Universe, M{pciit{z)), 
with the NFW density profile, expressed in units of enclosed 
mass and mean density within r, M{{p){< r)) at 2 = 0, in 
a mass-density plane. 

In this work we aim to provide a model that links the 
halo mass history with the halo concentration, a parameter 
that fully describes the internal structure of dark matter 
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halos. By doing so, we will gain insight into the origin of the 
NFW profile and its connection with the halo mass history. 
We also aim to find a physical explanation for the known 
correlation between the linear rms fluctnation of the density 
field, a, and halo concentration. 

This paper is organized as follows. We briefly introduce 
our simulations in Section[2l where we explain how we calcu¬ 
lated the merger history trees and discuss the necessary nu¬ 
merical convergence conditions. Then we provide a model for 
the halo mass history, which we refer to as the semi-analytic 
model. This semi-analytic model is described in Section [3l 
along with an analysis of the for mation time def i nition . For 
this model we use the empirical iMcBride et al.l l|200ilh for¬ 
mula. This function al form was motivate d by EPS theory in 
a companion paper (jCorrea et al ] l2015al: hereafter Paper I), 
and we calibrate the correlation between its two-parameters 
(a and /3) using numerical simulations. As a result, the 
semi-analytic model combines analytic relations with fits 
to simulations, to relate halo structure to the mass accre¬ 
tion history. In Section 13.61 we show how the semi-analytic 
model for the halo mass history depends on cosmology and 
the adopted definition for halo mass. In Section [4] we pro¬ 
vide a detailed comparison between the semi-analytic halo 
mass history model provided in this work, and the analytic 
model presented in Paper I. The parameters in this analytic 
model depend on the linear power spectrum and halo mass, 
whereas in the semi-analytic model the parameters depend 
on concentration and halo mass. We therefore combine the 
two models to establish the physical relation between the 
linear power spectrum and halo conce ntration. We will ex- 
pand on this in a forthcoming paper dCorrea et al.l[^15d . 
hereafter Paper III), where we predict the evolution of the 
concentration-mass relation and its dependence on cosmol¬ 
ogy. Finally, we provide a summary of formulae and discuss 
our main findings in Section [S] 


2 SIMULATIONS 

In this work we use the set of cosmological hydrodynam- 
ical simulations (the REF model) along with a set of 
dark ma tter only (DM ONLY) simulations from the OWLS 
project dSchave et al.l l2010lb These simulations were run 
with a significantly extended version of the N-Body Tree- 
PM, smoothed pa rticle hydrodyn amics (SPH) code Gadgets 
(last described in ISpringe]||2005l l. In order to assess the ef¬ 
fects of the finite resolution and box size on our results, 
most simulations were run using the same physical model 
(DMONLY or REF) but different box sizes (ranging from 
25/i“^Mpc to 400/i~^Mpc) and particle numbers (ranging 
from 128^ to 512®). The main numerical parameters of the 
runs are listed in Table 1. The simulation names contain 
strings of the form LxxxNyyy, where xxx is the simulation 
box size in comoving h“^Mpc, and yyy is the cube root 
of the number of particles per species (dark matter or bary- 
onic). For more details on the simulations we refer the reader 
to Appendix lAl 

Our DMONLY simulations assume the WMAP5 cos¬ 
mology, whereas the REF simulations assume WMAP3. To 
investigate the dependence on the adopted cosmological pa¬ 
rameters, we include an extra set of five dark matter only 
simulations (100 /i“®Mpc box size and 512® dark matter par- 


Table 2. Cosmological parameters. 


Simulation 

Om 

Ca 

h 

0-8 

Us 

DMONLY-WMAPl 

0.25 

0.75 

0.73 

0.90 

1.000 

DMONLY-WMAP3 

0.238 

0.762 

0.73 

0.74 

0.951 

DMONLY-WMAP5 

0.258 

0.742 

0.72 

0.796 

0.963 

DMONLY-WMAP9 

0.282 

0.718 

0.70 

0.817 

0.964 

DMONLY-Piancfcl 

0.317 

0.683 

0.67 

0.834 

0.962 


tides) which assume values for the cosmological parameters 
derived from the different releases of Wilkinson Microwave 
Anisotropy Probe (WMAP) and the Planck mission. Table 
[2] lists the sets of cosmological parameters adopted in the 
different simulations. 

Halo mass histories are obtained from the simulation 
outputs by building halo merger trees. We define the halo 
mass history as the mass of the most massive halo (main 
progenitor) along the main branch of the merger tree. The 
method used to create the merger trees is described in de¬ 
tail in Appendix \M\ While analysing the merger trees from 
the simulations, we look for a numerical resolution criterion 
under which mass accretion histories converge numerically. 
We begin by investigating the minimum number of particles 
halos must contain so that the merger trees lead to accurate 
numerical convergence. We find a necessary minimum limit 
of 300 dark matter particles, which corresponds to a mini¬ 
mum dark matter halo mass of Mhaio ~ 2.3 x 10®^ Mq in the 
100h“®Mpc box, Mhaio ~ 2.6 x 10®°Mq in the 50/i“®Mpc 
box, and Mhaio ~ 3.4 x 10®Mq in the 25/i“®Mpc box. 

In a merger tree, when a progenitor halo contains less 
than 300 dark matter particles, it is considered unresolved 
and discarded from the analysis. As a result, the number of 
halos in the sample that contribute to the median value of 
the mass history decreases with increasing redshift. Remov¬ 
ing unresolved halos from the merger tree can introduce a 
bias. When the number of halos that are discarded drops 
to more than 50% of the original sample, a spurious upturn 
in the median mass history occurs. To avoid this bias, the 
median mass history curve is only built out to the redshift 
at which less than 50% of the original number of halos con¬ 
tribute to the median mass value. 

Fig. [T] shows the effects of changing the resolution for 
the dark matter only and reference simulations. We first vary 
the box size while keeping the number of particles fixed (left 
panel). Then we vary the number of particles while keeping 
the box size fixed (right panel). The left panel (right panel) 
of Fig.[T]shows the mass history as a function of redshift for 
halos in eleven (seven) different mass bins for the DMONLY 
(REF) simulation. All halo masses are binned in equally 
spaced logarithmic bins of size A logj^Q M — 0.5. The mass 
histories are computed by calculating the median value of 
the halo masses from the merger tree at a given output red¬ 
shift, the error bars correspond to Icr confidence intervals. 
The different coloured lines indicate the different simulations 
from which the halo mass histories were calculated. The hor¬ 
izontal dash-dotted lines in the panels show the 300 x mam 
limit for the simulation that matches the colour. Halos in the 
simulation with masses lower than this value are unresolved, 
and hence their mass histories are not considered. The mass 
histories from halos whose main progenitors have masses 
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Table 1. List of simulations. From left-to-right the columns show: simulation identifier; comoving box size; number of dark matter 
particles (there are equally many baryonic particles); initial baryonic particle mass; dark matter particle mass; comoving (Plummer- 
equivalent) gravitational softening; maximum physical softening; final redshift. 


Simulation 

L 

(/i~^Mpc) 

N 

mb 

(h-iM0) 

^dm 

(/x-IMq) 

^com 

^prop 

(/i“^kpc) 


REF_L100N512 

100 

5123 

8.7 X 10’^ 

4.1 X 10® 

7.81 

2.00 

0 

REF_L100N256 

100 

2503 

0.9 X 10® 

3.2 X 10® 

15.02 

4.00 

0 

REF_L100N128 

100 

1283 

5.5 X 10® 

2.0 X lOi® 

31.25 

8.00 

0 

REF_L050N512 

50 

5123 

1.1 X 10’^ 

5.1 X 10'^ 

3.91 

1.00 

0 

REF_L025N512 

25 

5123 

1.4 X 10® 

0.3 X 10® 

1.95 

0.50 

2 

REF_L025N256 

25 

2503 

1.1 X 10’^ 

5.1 X 10'^ 

3.91 

1.00 

2 

REF_L025N128 

25 

1283 

8.7 X 10’^ 

4.1 X 10® 

7.81 

2.00 

2 

DMONLY_WMAP5_L400N512 

400 

5123 

- 

3.4 X lOi® 

31.25 

8.00 

0 

DMONLY-WMAP5_L200N512 

200 

5123 

— 

3.2 X 10® 

15.02 

4.00 

0 

DMONLY-WMAP5_L100N512 

100 

5123 

— 

5.3 X 10® 

7.81 

2.00 

0 

DMONLY-WMAP5-L050N512 

50 

5123 

— 

0.1 X 10'^ 

3.91 

1.00 

0 

DMONLY-WMAP5-L025N512 

25 

5123 

- 

8.3 X 10® 

2.00 

0.50 

0 



Log,o (1 +z) 



Log,o (1 +z) 


Figure 1. Median halo mass history as a function of redshift from simulations DMONLY (left panel) and REF (right panel) for halos in 
eleven and seven different mass bins, respectively. The curves show the median value, and the Icr error bars are determined by bootstrap 
resampling the halos from the merger tree at a given output redshift. The different colour lines show the mass histories of halos from 
different simulations. We find that a necessary condition for a halo to be defined, and mass histories to converge, is that halos should 
have a minimum of 300 dark matter particles. The horizontal dashed dotted lines show the 300 x mdm limit for the simulation that 
matches the colour, where mdm is the respective dark matter particle mass. When following a merger tree from a given halo sample, 
some halos are discarded when unresolved. This introduces a bias and so an upturn in the median mass history. Therefore, mass history 
curves are stopped once fewer than 50% of the original sample of halos are considered. Simulations in the REF model with 25/i“^Mpc 
comoving box size have a final redshift of z = 2, therefore the halos mass histories begin at this redshift. 
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lower than IO^^Mq at 2 = 0 were computed from simula¬ 
tions with 50/i“^Mpc and 25/i~^Mpc comoving box sizes. 
In the right panel, where the mass histories from the REF 
model are shown, all mass history curves obtained from the 
REF simulation with a 25h“^Mpc comoving box size have 
a final redshift of 2 = 2. Therefore, these halo mass histories 
begin at this redshift. 


3 SEMI-ANALYTIC MODEL FOR THE HALO 
MASS HISTORY 

In the following subsections we study dark matter halo prop¬ 
erties and provide a semi-analytic model that relates halo 
structure to the mass accretion history. We begin with the 
NEW density profile and derive an analytic expression for 
the mean inner density, (p)(< r- 2 ), within the scale radius, 
r_ 2 . We then define the formation redshift, and use the sim¬ 
ulations to find the relation between (p}(< r- 2 ) and the 
critical density of the universe at the formation redshift. We 
discuss the universality of the mass history curve and show 
how we can obtain an semi-analytic model for the mass his¬ 
tory that depends on only one parameter (as expected from 
our EPS analysis presented in a companion paper). We then 
calibrate this single parameter fit using our numerical sim- 


Table 3. Notation reference. Unless specified otherwise, quanti¬ 
ties are evaluated at 2 = 0. 


Notation 

Definition 

M 200 

Mr(r20o)-) total halo mass 

r200 

Virial radius 

r -2 

NEW scale radius 

c 

NFW concentration 

Mz 

M{z)^ total halo mass at redshift 2 : 

Mr{r) 

M(< r), mass enclosed within r 

X 

r/r-200 

(p)(< r- 2 ) 

Mean density within r _2 

Mr(r- 2 ) 

M(< r— 2 )i enclosed mass within r _2 

Z-2 

Formation redshift, when Mz equals Mr(r— 2 ) 

pcrit,0 

Critical density 

pcrit (-2) 

Critical density at redshift 2 : 

pm(^) 

Mean background density at redshift 2 : 

that the halo 

mass is defined as all matter within the radius 


^200 (see Table [ 3 ] for reference). 

The NEW profile is characterized by a logarithmic 
slope that steepens gradually from the centre outwards, 
and can be fully specified by the concentration param¬ 
eter and halo mass. Simulations have shown that these 
two parameters are correlated, with the average concentra- 


ulations. Finally, we show how the semi-analytic model for 

tion of a halo being a weakly decreasin 

g function of mass 

halo mass history depends on cosmology and halo mass def- 

fe.g. NFWilBullock et al.ll200ll: lEke et a 

.II 2 OOII: IShaw et al.l 

inition. 

200fil: Neto et al.l 20071: Duffv et al.l 20081: 

Maccio et al.ll2008l: 


Dutton & MaccidI I 2 OI 4 I: iDiemer & KravtsovI 1201,4'). There- 


3.1 Density profile 

An important property of a population of halos is their 
spheric ally averaged d e nsity profile. Based on Y-body simu¬ 
lations, I Navarro et al.l (|l997f l found that the density profiles 
of CDM halos can be approximated by a two parameter pro¬ 
file, 


p(r) = 


Pcrit^c 


{r/r- 2 ){l + r/r- 2 Y 


( 1 ) 


where r is the radius, r _2 is the characteristic radius 
at which the logarithmic density slope is —2, pcrit{z) = 
3H^{z)/8tvG is the critical density of the universe and Sc 
is a dimensionless parameter related to the concentration c 
by 


fore, the NEW density profile can be described by a single 
free parameter, the concentration, which can b e related to 
virial mass. The following relation was found bv IPuffv et al.l 
ll2008ll from a large set of Y—body simulations with the 
WMAP5 cosmology, 

c = 6.67(M2oo/2 X (3) 

for halos in equilibrium (relaxed). 

The NEW profile can be expressed in terms of the mean 
internal density 

Mr{r) _ 200 y(ca;) 


(p)(<r) = 


(47r/3)r 


Y{c) 


■Pcrit? 


(4) 


5c = 


200 


3 [ln(l + c) - c/(l -I- c)] ’ 


( 2 ) 


where x is defined as a: = r/r 2 oo and Y{u) = ln(l -|- u) — 
u/(l -I- u). Erom this last equation we can verify that at 
r = r2oo, X = 1 and (p)(< r 2 oo) = 200pcrit. 

Evaluating (p)(< r) at r = r_ 2 , we obtain 


which applies at fixed virial mass and where c is defined 
as c = r 2 oo/r- 2 , and r 2 oo is the virial radius. A halo is of¬ 
ten defined so that the mean density (p)(< r) within the 
halo virial radius ta is a factor A times the critical density 
of the universe at redshift 2 . Unfortunately, not all authors 
adopt the same definition, and readers should be aware of the 
difference in halo formation history and intern al structure 
when different mass defin itions are adopted Isee lDuffv et al.l 
l2008l : [D^iemer et al.ll2013l h We explore this in Section r3.6.2l to 
which the reader is referred to for further details. Through¬ 
out this work we use A = 200. We denote Mz = M 2 oo(z) as 
the halo mass as a function of redshift, Mr = M{< r) as the 
halo mass profile within radius r at 2 = 0, r 2 oo as the virial 
radius at 2 = 0 and c as the concentration at 2 = 0. Note 


(p)(< r_2) = 


Mr{r-2) 

A'^/i)C-2 


= 200c' 


1(1) 

Y(c) 


Pcrit • 


(5) 


From this last expression we see that for a fixed redshift the 
mean inner density (p}(< 'r- 2 ) can be written in terms of c. 
By substituting eq. into ®, we can obtain (p)(< r- 2 ) as 
a function of virial mass. Finally, we can compute the mass 
enclosed within r- 2 - From eq. m we obtain 

MrC- 2 ) = M 200 ■J', , , (6) 

Y(c) 

where we used M 200 = (47r/3)r2oo200pcrit. 

Although the NEW profile is widely used and gen¬ 
erally describes halo density profiles with high accuracy, 
it is worth noting that high resolution numerical simu¬ 
lations have shown that the spherically averaged density 
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L09io(^20o/^o) 


13 


14 

5 6 


15 

7 8 9 


DMONLY _ WMAP 5 _ L 400 N 51 2 [ Log , 0 ( 1 ^ 200 / 1 ^ 0 ) = 14 . 5-15 . 5 ] 

DMONLY_WMAP5_L2OON512[Log,o(M2oo/M0) = 13.5-14.5] 
DM0NLY_WMAP5_L100N51 2[Log,o(M2oo/Ms)= 1 3-14] 
DMONLY_WMAP5_L050N51 2[Log,o(M2oo/Ma) = 1 2-13] 
DMONLY_WMAP5_L025N512[Log,o(M2oo/Mg)=l l-i; 



■0> Relaxed sample - 
■0> Complete sample: 
<p>(<r_2) = 900porit(z_2) ; 
<p>(<r_2) = 854p„n(z_2) : 


0 . 0 - 0.5 1.0 1.5 

LOgio Pcrlt(Z-2)/Pcrlt,0 


2.0 


Figure 2. Relation between the mean density within the NFW 
scale radius at 2 : = 0 and the critical density of the universe at the 
halo formation redshift, Z— 2 -, for DMONLY simulations from the 
OWLS project. The simulations assume the WMAP-5 cosmolog¬ 
ical parameters and have box sizes of 400h“^Mpc, 200/i“^Mpc, 
100h“^Mpc, 50/i“^Mpc and 25/i“^Mpc. The black solid line in¬ 
dicates the relation shown in eq. which only depends on cos¬ 
mology through the mass-concentration relation. The black (red) 
star symbols show the mean values of the complete (relaxed) sam¬ 
ple in logarithmic mass bins of width 51ogjQ M = 0.2. The black 
dashed and solid lines show the relations found by fitting the data 
of the complete and relaxed samples, respectively. The filled cir¬ 
cles correspond to values of individual halos and are coloured by 
mass according to the colour bar at the top of the plot. 


profiles of dark matter halos have s mall but systematic 
deviations from the N F W form (e.g. [Navarro et al. 200j; 
Hwashi&^^hi^ 200^ Narorro et ahlhoid : iLudlow et al.l 


2OI0I : Diemer fc Kravtsov! 20 l 4 ] . While there is no clear un¬ 


derstanding of what breaks the structural similarity among 
halos, an alternative parametrization is sometimes used (the 
Einasto profile), which assumes the logarithmic slop e to be a 
simp le power la w of radius, dIn p/d In r oc (r/r_ 2 )“ llEinastd 
Il96d) . Recently. iLudlow et al.l (l2013l ) investigated the rela¬ 
tion between the accretion history and mass profile of cold 
dark matter halos. They found that halos whose mass pro¬ 
files deviate from NFW and are better approximated by 
Einasto profiles also have accretion histories that deviate 
from the NFW shape in a similar way. Ffowever, they found 
the residuals from the systematic deviations from the NFW 
shape to be smaller than 10%. We therefore only consider 
the NFW halo density profile in this work. 


3.2 Formation redshift 

iNavarro et al.l (Il997l ) showed that the characteristic over¬ 
density (5c) is closely related to formation time (zf), which 
they defined as the time when half the mass of the main pro¬ 
genitor was first contained in progenitors larger than some 
fraction / of the mass of the halo at 2 = 0. They found 
that the ‘natural’ relation 5c oc flm(l + zt)^ describes how 


the overdensity of halos varies with their formation redshift. 
Subsequent investigations have used N-body simulations and 
empirical models to explore the relation bet ween concentra¬ 
tion and formation history i n more detail llWechsler et al.l 
I2OO2I: IZhao et aklbool I2009I) . A good definition of forma¬ 
tion time that relates concentration to halo mass history 
was found to be the time when the main progenitor switches 
from a period of fast growth to one of slow growth. This is 
based on the observation that halos that have experienced a 
recent major merger typically have relatively low concentra¬ 
tions, while halos that have experienced a longer phase of rel- 
ative ly quiescent growt h have larger concentrations. More¬ 
over, (zhao et ^ (I2OO9I) argue that halo concentration can 
be very well determined at the time the main progenitor of 
the halo has 4% of its final mass. 

The various formation time definitions each provide ac¬ 
curate fits to the simulations on which they are based and, 
at a given halo mass, show reasonably small scatter. How¬ 
ever, our goal is to adopt a formation time definition that 
has a natural justification without invoking arbitrary mass 
fractions. To this end, we go back to the idea that halos are 
formed ‘inside out’, and consider the formation time to be 
defined as the time when the initial bound core forms. We 
follow Ludlow et al. (2013) and define the formation redshift 
as the time at which the mass of the main progenitor equals 
the mass enclosed within the scale radius at 2 = 0, yielding 

2:_2 = z[M^ = Mr{r- 2 )\. (7) 

From no w on we denote the fo rmation redshift by Z- 2 - Inter¬ 
estingly, ILudlow et al.l (I 2 OI 3 I ) found that at this formation 
redshift, the critical density of the universe is directly pro¬ 
portional to the mean density within the scale radius of halos 
at 2 = 0 : pcrit(z- 2 ) oc (p)(< r- 2 ). A possible interpretation 
of this relation is that the central structure of a dark matter 
halo (contained within r- 2 ) is established through collapse 
and later accretion and mergers increase the mass and size of 
the halo without adding much material to its inner regions, 
thus increasing the halo virial radius while leaving the scale 
radius and its inner density {(p){< r- 2 )) almost unchanged 
||Huss et al.lll999l : IWang fc Whitdl2009l ). 


3.3 Relation between halo formation time and 
concentration from numerical simulations 

We now study the relation between pcrit( 2 - 2 ) and (p)(< r- 2 ) 
using a set of DMONLY cosmological simulations from the 
OWLS project that adopt the WMAP-5 cosmology. We be¬ 
gin by considering two samples of halos. Our complete sam¬ 
ple contains all halos that satisfy our resolution criteria while 
our ‘relaxed’ sample retains only those halos for which the 
separation between the most bound particle and the centre 
of mass of the Friends-of-Friends (FoF) halo is smaller than 
0.077?vir, where fJvir is the ra dius for which the mean inter¬ 
nal density is A (as given by Brv an fc N ormarjn~998l) times 
the critical density. iNeto et al.l ( 20071 ) found that this sim¬ 
ple criterion resulted in the removal of the vast majority of 
unrelaxed halos and as such we do not use their additional 
criteria. At z = 0, our complete sample contains 2831 halos, 
while our relaxed sample is reduced to 2387 (84%). 

To compute the mean inner density within the scale ra¬ 
dius, (p)(< r- 2 ), we need to fit the NFW density profile 
to each individual halo. We begin by fitting NFW profiles 
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.I.I.11111 

ODMONLY.WMAP5_L400N512 ^ 
ADMONLY.WMAP5_L200N512 ! 



- - . (1 +z.2)^=[(200/854)c^ Y(1 )/Y(cc)-Q v 
(1+z.,)^=[(200/900)c^Y(1)/Y(c„)-QjQ-’ ; 


.I.I.I.I I I I I 

11 12 13 14 15 

Log^o Mq [Mq] 



c 


Figure 3. Relation between formation redshift, Z— 2 , and halo concentration, c (left panel), and between formation redshift and z = 0 
halo mass, M 200 (right panel). The different symbols correspond to the median values of the relaxed sample and the error bars to Icr 
confidence limits. The solid line in the left panel is not a fit but a prediction of the Z —2 — c relation for relaxed halos given by eq. (El). 
Similarly, the dashed line is the prediction for the complete sample of halos, assuming a constant of proportionality between (p)(< r— 2 ) 
and pcrit of 854, rather than the value of 900 used for the relaxed sample. The grey area shows the scatter in z —2 plotted in Fig. lAll fright 
panel). Similarly, the solid line in the right panel is a prediction of the z —2 — M 200 relation given by eqs. m and 0- The dashed line 
also shows the z —2 — M 200 relation assuming (p)(< r— 2 ) = 854perit a-nd the concentration-mass relation calculated using the complete 
sample. 


to all halos at 2 : = 0 that contain at least 10^ dark matter 
particles within the virial radius. For each halo, all particles 
in the range —1.25 ^ log,Q(r/r 2 oo) ^ 0, where r 2 oo is the 
virial radius, are binned radially in equally spaced logarith¬ 
mic bins of size A log,g r = 0.078. The density profile is then 
fit to these bins by performing a least square minimization 
of the difference between the logarithmic densities of the 
model and the data, using equal weighting. The correspond¬ 
ing mean enclosed mass, Mr(r_ 2 ), and mean inner density 
at r_ 2 , (p)(< r- 2 ), are found by interpolating along the cu¬ 
mulative mass and density profiles (measured while fitting 
the NFW profile) from r = 0 to r _2 = rzoo/c, where c is the 
concentration from the NFW fit. Then we follow the mass 
history of these halos through the snapshots, and interpolate 
to determine the redshift Z -2 at which = Mr{r- 2 ). 

We perform a least-square minimization of the quan¬ 
tity A^ = ^ y;,^,[{pi)(Pcrit,i) - /(pcrit,i,4l)], to obtain 

the constant of proportionality, A. We find (p)(< r- 2 ) = 
(900±50)pcrit(2-2) for the relaxed sample, and (p)(< r_ 2 ) = 
(854 ± 47)pcrit(2-2) for the complete sample. The Icr er- 
ror was obtained fro m the least squares fit. For comparison, 
iLudlow et al.l (l2014l ~l found a constant value of 853 for their 
relaxed sample of halos and the WMAP-1 cosmology. Fig. [2] 
shows the relation between the mean inner density at 2 = 0 
and the critical density of the universe at redshift 2_2 for 
various DMONLY simulations. Each dot in this panel cor¬ 
responds to an individual halo from the complete sample 
in the DMONLY-WMAP5 simulations that have box sizes 
of 400/i~^Mpc, 200/i“^Mpc, 100h“^Mpc, 50h“^Mpc and 


25/i“^Mpc. The (p)(< r- 2 )-pcrit{z- 2 ) values are coloured 
by mass according to the colour bar at the top of the plot. 
The black (red) star symbols show the mean values of the 
complete (relaxed) sample in logarithmic mass bins of width 
51og,o M = 0.2. As expected when unrelaxed halos are dis¬ 
carded (e.g. Duffy et al. 2008), the relaxed sample contains 
on average slightly higher concentrations (by a factor of 1.16) 
and so higher formation times (by a factor of 1.1). 

In Fig. [2]the best-fit to the data points from the relaxed 
sample is shown by the solid line, while the dashed line is 
the fit to the complete sample. The Pcrit(z- 2 ) — (p)(< r- 2 ) 
correlation clearly shows that halos that collapsed earlier 
have denser cores at 2 = 0. Using the mean inner density- 
critical density relation for the relaxed sample, 

(p)(<r_ 2 ) = (900±50)pcrit(2-2). (8) 


We replace (p)(< r- 2 ) by eq. ([5]) and calculate the depen¬ 
dence of formation redshift on concentration. 


200 Y(l) Da 


(9) 


This last expression is tested in Fig. [3] (left panel) where we 
plot the median formation redshift as a function of concen¬ 
tration using different symbols for different sets of simula¬ 
tions from the OWLS project. The symbols correspond to 
the median values of the relaxed sample, and the error bars 
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indicate la confidence limits. The grey solid line shows the 
Z-2 — C relation given by eq. ( 0 ) , whereas the grey dashed line 
shows the same relation assuming a constant of 854 instead 
of 900 (as obtained for the complete sample). Similarly, us¬ 
ing the Duffy et al. (2008) concentration-mass relation we 
obtain the formation redshift as a function of halo mass at 
z = 0 (right panel of Fig. [3]). It is important to note that the 
Z -2 —c and Z -2 — M relations are valid in the halo mass range 
10 ^^ — 10^® at lower masses the conce ntration-mass 

relati on begins to deviate from power laws (iLudlow et al.l 
l2014h . 

In Appendix we analyse the scatter in the forma¬ 
tion time—mass relation and show that it correlates with 
the scatter in the concentration—mass relation. Thus con¬ 
cluding that the scatter in formation time determines the 
scatter in the concentration. Also, we investigate how the 
scatter in halo mass history drives the scatter in formation 
time. 


3.4 The mass history 

Fig. 2] shows the mass accretion history of halos in differ¬ 
ent mass bins as a function of the mean background density. 
The mass histories are scaled to Mr{r- 2 ) and the mean back¬ 
ground densities are scaled to pm(z- 2 ) = Pcrit.ol^ m (1 + ^-2)®. 
The figure shows that all halo mass histori es look alike. This 
is in agreement with iLudlow et al.l (l2013f l , who found that 
the mass accretion history, expressed in terms of the critical 
density of the Universe, M(pcrit(z)), resembles that of the 
enclosed NFW density profile, M((p)(< r)). The similarity 
in the shapes between M{pcrit{z)) and M{{p){< r)) is still 
not clear, but it suggests that the physically motivated form 
M{z) = Mo(l -I- z)°‘e^^, which is a result of rapid growth in 
the matter dominated epoch followed by a slow growth in 
the dark energy epoc h, produces the double power-law of the 
NFW profile (see e.g. ILu et al.ll200(^ . We use this feature to 
find a functional form that describes this unique universal 
curve in order to obtain an empirical expression for the mass 
accretion history at all redshifts and halo masses. 

We are motivated by the extended Press-Schechter anal¬ 
ysis of halo mass histories presented in Paper I. In that work, 
we showed through analytic calculations, that when halo 
mass histories are described by a power-law times an expo¬ 
nential, 

M(z) = Mo(l-f z)“e^", (10) 

and that the parameters a and /I are connected via the power 
spectrum of density fluctuations. In this work however, we 
aim to relate halo structure to the mass accretion history. 
We therefore first determine the correlation between the pa¬ 
rameters a and /3 and concentration. To this end, we first 
find the a — j3 relation that results from the formation red- 
shift definition discussed in the previous section. Thus, we 
evaluate the halo mass at Z-2, 


M.{z- 2 ) = Mr{r. 2 ) = M.(« = 0 )(l-bz- 2 )“e^"-A (11) 


Taking the natural logarithm, we obtain, 

= ^H^ + Z-2) + Pz.2, (12) 



In Pm(z)/Pm(Z-2) 

Figure 4. Mass histories of halos, obtained from different 
DMONLY_WMAP5 simulations, as indicated by the colours. 
The bottom left legends indicate the halo mass range at 2 = 0, 
selected from each simulation. For example, we selected halos be¬ 
tween 10® — from the DMONLY_WMAP5_L025N512 

simulation, divided them in equally spaced logarithmic bins of 
size Alogj^Q M = 0.2, and calculated the median mass histories. 
The different curves show the median mass history of the main 
progenitors, normalized to the median enclosed mass of the main 
progenitors at 2 : = 0, My.(r_ 2 ). The mass histories are plotted as 
a function of the mean background density of the universe, scaled 
to the mean background density at Z— 2 - The blue dashed line is 
a fit of expression in to the different mass history curves. The 
median value of the only adjustable parameter, 7 , is indicated in 
the top-right part of the plot. 


and hence 



ln(l -I- Z-2) 


(13) 


From this last equation we see that a can be written as a 
function of /3, Mr{r- 2 ), M^{z = 0) and Z- 2 - However, as 
Mr{r- 2 ) is a function of concentration and virial mass (see 
eq.|6j, we can write a in terms of P, concentration and Z- 2 , 


^ HY{l)/Y{c))-pz -2 

\n{l + z- 2 ) 

The next step is to find an expression for p. We find 
P{z- 2 ) by fitting eq. (IIOII to all the data points plotted in 
Fig.|4] We now need to express M{z) leg. HOI) as a function of 
the mean background density. To do this, we replace (1 -|- z) 
by (pm(2)/pcrit,o/I2m)'^^® and divide both sides of eq. (fTOt 
by Mr{r- 2 ), yielding 


M,{z) 

Mr{r-2) 


M^{Z = 0) / Pm(z) 
Y[r{T — 2) yUniPcrit.O J 


X 


exp P 


pmjz) 

IlmPcrit.O 



(15) 


Multiplying both denominators and numerators by Pmiz- 2 ), 
we get, after rearranging. 
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MRz) 

Mr{r-2) 


Mz{z = 0) /f pmjz) 

— A Y^mPcritjO J yPm(^—2) J 


X exp P 


Pin(g-2) 

^mPcrit,0 


1/3 


- 1 


X 


exp 7 


f \ \ 
Pm(^) \ 

Pm(«-2) J 



(16) 


where, we have 

/3(pm(^-2)/flm/Pcrit.o)^^^ = 


defined 7 
P{1 + 2 - 2 ). The 


term 


^ Pm(z- 2 ) \ 
y flmPcritjO J 


M^(z=0) _ 

Mr(r_2) \ nmPorit,( 


a/3 


exp 


f Pin(z- 2 ) \ 
r2niPcrit,0 / 


1/3 


- 1 


in eq. (nsj is equal to unity, which can be seen by replacing 
pui[z- 2 )/ 0 .rapciitfi = (1 + 2 - 2)3 and comparing with eq. 
dj. Hence eq. (Hi) becomes 


MRz) 

Mr{r-2) 


X 


Pm(2) 

Pm(2_2) 

exp I 7 


y./3 


Pm{z) 

PmiZ- 2 ) 



(17) 


Thus, based on eq. mi), the functional form to fit the 
mass accretion histories from the simulations can be written 
as 


/( 2 , 7 ) = a( 2 - 2 , c, 7 ) 2/3 + 7 (e^/® - 1 ), (18) 


where 7 ( 2 , 7 ) = In (and 2 = In 

/ \Mr(r_ 2 ) J [Pm(2-2 


. From 


eq. dll we see that the parameter a is a function of z- 2 , c 
and 7 , 


^ ^ ln(y(l)/y(c)) - 72 - 2/(1 + 2 - 2 ) 

ln(l + 2 _ 2 ) 

Therefore, 7 is now the only adjustable parameter. We per¬ 
form a x^-like minimization of the quantity 


N 

^ '^[logio{Mz{zi)/Mrir-2)) - /(2i,7)l^ (20) 

i=l 


and find the value of 7 that best fits all halo accretion his¬ 
tories. The sum in the x^-like minimization is over the N 
available simulation output redshifts at Zi{i = 1,N), with 

Zi = In 


Pcrit,0^m(l-t-2i)'^ 


Pm(2-2) 


Fig- HI shows halo mass histories (with M,,{z) scaled to 
Mr{r- 2 )) for our complete halo sample as a function of the 
mean background density [pm( 2 ) scaled to pm( 2 - 2 )]. The 
blue dashed line is the fit of expression m to all the mass 
history curves included in the figure. Here, the only ad¬ 
justable parameter is 7 . We obtained 7 = —3.01 ± 0.08, 
yielding 


P = - 3 (p^( 2 - 2 )/fl,,/Pcrit, 0 )"^^® = -3/(1 -b 2 - 2 ). (21) 


Fig. |4] shows that the halo mass histories have a char¬ 
acteristic shape consisting of a rapid growth at early times. 


followed by a slower growth at late times. The change from 
rapid to slow accretion corresponds to the transition between 
the mass and dark energy dominated eras (see Paper I), and 
depends on the parameter P in the exponential (as can be 
seen from eq. HQl). The dependence of P on the formation 
redshift is given in eq. m, which shows that a more recent 
formation time, and hence a larger halo mass, results in a 
larger value of p, and so a steeper halo mass history curve. 
This last point can be seen in Fig.[5]from the mass histories 
of halos in different mass bins (coloured lines shown in the 
panels). The panel on the left shows the mass history curves 
from the DMONLY simulation outputs (coloured lines in the 
background), and the mass histories predicted by eqs. (HOI) , 
m and (I 2 III (red dashed lines). From these panels we see 
that, (i) the mass history formula works remarkably well 
when compared with the simulation, and (ii) the larger the 
mass of a halo at 2 = 0 , the steeper the mass history curve at 
early times. In contrast, the mass history of low-mass halos 
is essentially governed by the power law at late times. 

The halo mass histories plotted in Fig.[5]come from the 
complete sample of halos (relaxed and unrelaxed). We found 
no significant difference in mass growth when only relaxed 
halos are considered. We therefore conclude that the fact 
that a halo is unrelaxed at a particular redshift does not af¬ 
fect its halo mass history, provided the concentration—mass 
relation fit from the relaxed halo sample is used. This is an 
interesting result because while deriving the semi-analytic 
model of halo mass history, we assumed that the halo density 
profile is described by the NFW profile at all times. There¬ 
fore while the N FW is not a good fit for the density profile of 
unrelaxed halos (iNeto et al.ll2007ll . our semi-analytic model 
(based on NFW profiles) is a good fit all halos (relaxed and 
unrelaxed). 

The right panel of Fig.[5]shows halo mass histories from 
the REF hydrodynamical simulations. We compute the halo 
mass as the total mass (gas and dark matter) within the 
virial radius (r 2 oo). We find that the inclusion of baryons 
steepens the mass histories at high redshift, therefore the 
best description of M{z) is given by eqs. HOII . H4I) . (I21II . and 
the concentration-mass relation from the complete sample of 
halos, c = 5.74(M/2 x 


3.5 The mass accretion rate 

The accretion of gas and dark matter from the intergalac- 
tic medium is a fundamental driver of both, the evolution of 
dark matter halos and the formation of galaxies within them. 
For that reason, developing a theoretical model for the mass 
accretion rate is the basis for analytic and semi-analytic 
models that study galaxy formation and evolution. In this 
section we look for a suitable expression for the mean accre¬ 
tion rate of dark matter halos. To achieve this, we take the 
derivative of the semi-analytic mass history model, M{z), 
given by eq. (HSl) with respect to time and replace dz/dt by 
-Ho[nin(l + zf + flA(l + yielding 

= 71.6MQyr“^Mi2/io.7[-a -/3(1-I- 2 )] 

X [IIm(l + 2)3-1 -Ha] ^ , (22) 

where ho .7 = hjQ.l, M 12 — M/IO^^Mq and a and P are 
given by eqs. and dnj, respectively. As shown in the 
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Log,o(1+z) Log,o(1+z) 


Figure 5. Mass histories of all halos from simulations DMONLY_WMAP5 (left panel) and REF (right panel). Halo masses are binned 
in equally spaced logarithmic bins of size A logj^Q M = 0.5. The mass histories are computed by calculating the median value of the halo 
masses from the merger tree at a given output redshift, the error bars correspond to Icr. The different colour lines show the mass histories 
of halos from different simulations as indicated in the legends, while the red dashed curves correspond to the mass histories predicted by 
eqs. fTOll. l|T9ll and (IHJ. 


previous section, the parameters a and /3 depend on halo 
mass (through the formation time dependence). We hnd that 
this mass dependence is crucial for obtaining an accurate de¬ 
scription for the mass history (as shown in Fig. However, 
the factor of 2 (3) change for a (/3) between halo masses 
of 10® and 10^^ Mq is not significant when calculating the 
accretion rate. Therefore, we provide an approximation for 
the mean mass accretion rate as a function of redshift and 
halo mass, by averaging a and /3 over halo mass, yielding 
{q)= 0.24, {/3H-0.75, and 

i-^) = 71.6MQyr“^Mi2/io.7 

X[-0.24 + 0.75(1 -b 2)][flm(l + zf + flA]^/f23) 

Fig. (top panel) compares the median dark matter 
accretion rate for different halo masses as a function of red¬ 
shift (solid lines) to the mean accretion rate given by eq. (1231) 
(grey dashed lines). From the merger trees of the main halos, 
we compute the mass growth rate of a halo of a given mass. 
We do this by following the main branch of the tree and 
computing dM/dt = (M(zi) — M(z 2 ))/At, where zi < Z 2 , 
M{zi) is the descendant halo mass at time t and M{z 2 ) is 
the most massive progenitor at time t — At. The median 
value of dM/dt for the complete set of resolved halos is then 
plotted for different constant halo masses. We find very good 
agreement between the simulation outputs and the analytic 


estimate given by eq. (1231) . As expected, the larger the halo 
mass, the larger the dark matter accretion rate. 

3.5.1 Baryonic accretion 

Next, we estimate the gas accretion rate and compare 
our model wi t h sim ilar fitting formulae propose d by 
iFakhouri et al.l (l2010ll and IPekel fc Krumhola (l2013l l . The 
bottom panel of Fig. [ 5 ] shows the gas accretion rate as a 
function redshift for a range of halo masses (log,g M /Mq = 
11.2 — 12.8). The grey circles correspond to the gas ac¬ 
cretion rate measured in REF_L100N512. In this case we 
compute the total mass growth {M — Mgas -I- Mum) from 
the merger trees, and then estimate the gas accretion rate 
by multiplying the total accretion rate by the universal 
baryon fraction fi, — flb/Pm. The green solid line corre¬ 
sponds to our gas accretion rate model (given by flb/Pm 
times eq. US]). The b lue dot-dashed line is th e gas accre¬ 
tion rate proposed bv IPekel fc Krumholj (l2013tl {dMh/dt = 
fbMi 2 {l -|- 2 )®^^), who derived the baryonic in¬ 
flow on to a halo dMh/dt from the averaged growth rate 
of halo mass through mergers and smooth a ccretion based 
on the EPS theory of gravita tional clustering (iNeistein et al.l 
I 2 OO 6 I : INeistein fc Pek^l2008ll . Lastly, we compare our mode l 
with the accretion rate formula from IFakhouri et al.l (I 2 OI 0 II 
[dMh/dt = 46.1M0yr“VbA7iV(l + l.ll2)(Pm(l + z)^ -|- 
Pa)^^ ^L plotted as the purple dashed line. IFakhouri et al.l 
1 I 2 OIOI) constructed merger trees of dark matter halos and 
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Log,o(1 +z) 



Log,o(1 +z) 


Figure 6. Mean accretion rate of dark matter (top panel) and 
gas (bottom panel) as a function of redshift for different halo 
masses. Top panel: the accretion rate obtained from the simu¬ 
lation outputs up to the redshift where the halo mass histories 
are converged. Grey dashed lines show the accretion rate esti¬ 
mated using eq. ll23t . Bottom panel: gas accretion rate obtained 
from the REF_L100N512 simulation (grey circles), from Slb/nin 
times eq. 123} (green solid line), and from various fitting formulae 
taken from the literature. 


quantified their merger rates and mass growth rate using 
the Millennium and Millennium II simulations. They de¬ 
fined the halo mass as the sum of the masses of all subhalos 
within a FoF halo. We see that our accretion rate model is in 
excell ent ag reement with the fo r mulae from iFakhouri et al.l 
1 I 2 OIOI) and Dekel fc K rumholj (l2013l ). We find that the 
IFakhouri et al.l ( 2010l) formula generally overpredicts the gas 
accretion rate in the low-redshift regime (e.g. it overpredicts 


it b y a factor of 1.4 at 2 = 0 for a IO^^Mq mass halo). 
The lPekel fc Krumholj l2013l) formula underpredicts (over¬ 


predicts) the gas accretion rate in the low- (high-) redshift 
regime for halos with masses larger (lower) than IO'^^Mq. 


3.6 Dependence on cosmology and mass definition 

We have developed a semi-analytic model that relates the 
inner structure of a halo at redshift zero to its mass history. 
The model adopts the NFW profile, computes the mean in¬ 
ner density within the scale radius, and relates this to the 
critical density of the universe at the redshift where the halo 
virial mass equals the mass enclosed within r_ 2 . This re¬ 
lation enables us to find the formation redshift-halo mass 
dependence and to derive a one-parameter model for the 
halo mass history. In this section we consider the effects of 
cosmology and mass definition on the semi-analytic model. 


3.6.1 Cosmology dependence 


The adopted cosmological parameters affect the mean in¬ 
ner halo densities, concentrations, formation redshifts and 
halo mass histories. To investigate the dependence of halo 
mass histories on cosmology, we have run a set of dark mat¬ 
ter only simulations with different cosmologies. Table[2]lists 
the sets of cosmological parameters adopted by the different 
simulations. Specifically, we assume values for the cosmo¬ 
logical parameters derived from measurements of the cos¬ 
mic micro wave background by the WMAP and the Planck 
missions j[Spergel et alJj200^!j2007];jKomnt^^ 


iHinshaw et al. 20131 : Planck Collaboration et al. 20141 ) . 

It has been sh own that halos that formed earlie r are 
more concentrat ed jNayarro et ah 1997|; BuIIoc^_e£_alJ 


Eke et al.l 


Neto et ^ 


2001 


lKuhIeii^£aI 


2007 ). iMaccio et al. 


2005 


iMaccib et al. 


2001 


2007 


( 20081 ) explored the depen¬ 


dence of halo concentration on the adopted cosmological 
model for field galaxies. They found that dwarf-scale field 
halos are more concentrated by a factor of 1.55 in WMAPl 
compared to WMAP3, and by a factor of 1.29 for cluster¬ 
sized halos. This reflects the fact that halos of a fixed 2 = 0 
mass assemble earlier in a universe with higher flm, higher 
as and/or higher Us. 

The halo formation redshift can be related to the power 
at the corresponding mass scale, and therefore depends on 
both (7g and n^. The parameter as sets the power at a scale 
of 8/i“^Mpc, which corresponds to a mass of about 1.53 x 
10^^/i“^MQ(IIin/fIm,WMAP5), and a waveuumber of fcg- This 
last quantity is given by the relation M = (47rpm/3)(27r/fc)®. 
For a power-law spectrum P(k) oc fc", the variance can be 
written as a^{k)/a^ = (fc/fcg)"^®. Therefore, the change in 
a between WMAP5 and WMAPl for a given halo mass that 
corresponds to a wavenumber fc is 


CrWMAPl(fc) ^ as,WMAPl ( fc WMAPl-» s,WMAP5)/2 

awMAP5(fc) ag.WMAPS 


(s) 


.(24) 


A halo mass of IO^^Mq corresponds to a wavenumber of 
fci .3 ~ Bfcg. The total change in the mean power spectrum at 
this mass scale is °'wMAPi(fei.3) _ x 27. This is proportional 

CTWMAPSI^I.S) 

to the change of the formation redshift, 

(1 + 2:f,WMAPl) = 1.27(1 -I- 2f,WMAP5)- (25) 

Next, we test how this change affects the halo mass 
history. We showed that the mass history profile is well 
described by the expression M{z) = Mz{z = 0)(1 + 2 )“e^^, 
where a and p both depend on the formation redshift. In 
the mass history model presented in Section 4 there are two 
best-fitting parameters that can be cosmology dependent. 
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One is the constant value A = 900 in eq. © that relates the 
mean inner density to the critical density at Z- 2 , and the 
other is the constant value 7 = —3 in eq. Isa that defines 
the P parameter. 

To investigate the cosmology dependence of j 4, we anal¬ 
yse the (p)(< r- 2 ) — pcrit{z- 2 ) relation in the simulations 
with different cosmologies. We do the same as in Section 4.3. 
First we fit the NFW profile to dark matter halos to obtain c 
and r- 2 , and calculate the cumulative mass, M- 2 , and den¬ 
sity, (p)(< r- 2 ), from r = 0 to r = r- 2 . Then we follow the 
halo mass histories through the snapshots and interpolate to 
calculate Z- 2 , the redshift for which M{z) is equal to M- 2 . 
Finally, we obtain the best-fitting {p){< r- 2 ) — pcrit{z- 2 ) 
relation. We find that the parameter Acosmo, where cosmo is 
WMAPl, WMAP3, WMAP5, WMAP9 or Planck, changes 
with cosmology. We show this in the top panel of Fig. [ 8 ] 
We find Awmapi = 787 ± 52.25, Awmaps = 850 ± 39.60, 
Awmaps = 903 ± 48.63, Awmapq = 820 ± 51.03 and 
Apia nck = 798 ± 43.73. We do not find good agreement 
with iLudlow et ^ (l2013l i. who found Awmapi = 853 for 
WMAPl cosmology. This is due to the fact that we are only 
analysing the {p){< r- 2 ) — pcrit(z- 2 ) relation in the high- 
mass regime (M = 10^^ ® — 10^^ ®Mq), due to the limitations 
of the box size (L = 100/i“'^Mpc). With a more complete 
halo population, we may obtain better agreement. We con¬ 
clude that the Acosmo parameter depends on cosmology, at 
least for the halo mass ranges we are considering. 

Next, we analyse how the change of the formation red- 
shift due to cosmology affects /3, defined as /3 = —3/(1 -I- z{). 
We find from Fig. [7] that the constant value, —3, is insensi¬ 
tive to cosmology. Fig. [3 shows the same analysis as Fig. [I] 
but for halo mass histories obtained from simulations with 
the WMAPl and WMAP5 cosmology as indicated in the leg¬ 
ends. We fit expression m to the mass history curves from 
different cosmologies and obtained the same adjustable pa¬ 
rameter 7 . We therefore conclude that 7 = —3 is insensitive 
to cosmology. 

If we then consider that the change in /? between 
WMAPl and WMAP5 is Pwmapi = —3/(1 -I- Z{, , WMAPl) — 
-3/1.27/(l-|-2f,wMAP5) = PwMAPs/i.27, the change in the 
halo mass history between the WMAP5 and WMAPl cos¬ 
mologies, for a halo mass of 10 ^^Mq at 2 = 0 , corresponds 
to 


Irio- ^(^)wMAPl _ , (/SwMAPl—PwMAPsl^i 

^1°M(z)wMAP5 

« —0.12/3wmapi2, 

« O.lz. (26) 


In the last step we replaced /3wmapi by 
—3/(1-I-WMAPl) =—0.75 for a IO'^^Mq halo. We 
obtained Zf,wMAPi from the {p){< r- 2 ) — pcAt{z-2) rela¬ 
tion suitable for the WMAPl cosmology (see t op panel of 
Fig. [3 ) and the concentration-mass relation from iNeto et al.l 

(l2003). 


Next, we test the change in halo mass history. For ex¬ 
ample, if a halo had a mass of 10 ^^ ‘^Mq at 2 = 2 in the 
WMAP5 cosmology, it would have had a mass of 10^^ ®Mq 
in the WMAPl cosmology. The value of as has a particu¬ 
larly large effect at high redshift, because structure forma¬ 
tion proceeds faster in the WMAPl cosmology, as shown by 
the above expression. This last point can also be seen in the 



In Pm(z)/Pm(Z-2) 

Figure 7. Mass histories of halos, obtained from simulations 
with the WMAPl cosmology (DMONLY_WMAP1_L100N512, 
blue solid lines) and the WMAP5 cosmology 
(DMONLY_WMAP5_L100N512, green solid lines). The 
curves show the median mass history of the main progenitors, 
normalized to the median enclosed mass, Mr(‘r_ 2 ), of the main 
progenitors at 2 = 0. The mass histories are plotted as a function 
of the mean background density of the universe, scaled to the 
mean background density at 2 _ 2 . The blue dashed line is a 
fit of expression GHJ to the different mass history curves. The 
median value of the only adjustable parameter, 7 , is indicated 
in the top-right part of the plot. We find that 7 is insensitive to 
cosmology. 


two panels of Fig. [S] The top panel is the same as the right 
panel of Fig. [31 and shows the formation redshift, Z- 2 , as a 
function of halo mass (obtained from simulations with dif¬ 
ferent cosmologies). We see that there are large differences 
between the WMAP5 and WMAPl cosmologies due to the 
changes in crs and ris. Interestingly, there is only a small dif¬ 
ference between the Planck and W MA Pl cosmologies (in 
agree ment with [Ludlow et al.l |2014| and [Dutton fc Macci^ 
l2014h . and also between the Planck and WMAP9 cosmolo¬ 
gies for which we found (l-|-2f)pianck = 1.01(1-1- 2f)wMAP9 for 
a IO^^Mq halo. The bottom panel of Fig. [Hjshows the mass 
history of 10^^ Mq halos at 2 = 0 from DMONLY simula¬ 
tions with different cosmologies. As predicted, the change 
in mass history between the WMAP5 and WMAPl cos¬ 
mologies is log;,Q M(z)wMAP 5 O.I 2 , while little difference is 
found between the WMAP9 and Planck {AM{z) ^ 10~^z). 

We found that as long as a suitable concentration-mass 
relation and value for the A parameter are assumed for the 
cosmology being considered, eqs. (l37ll . (i38l) and (l39l) pro¬ 
vide a good estimate of the mass history curve. This can be 
seen in the bottom panel of Fig. [8]by comparing the different 
halo mass histories. For the WMAP5 and WMAP3 cosmolo- 
gies we assumed t he concentration-mass relation found by 
jPuffy et al.l (l2008l ). whereas we used the relation from Neto 
et al. (2007) for the other cosmologies. For a step-by-step 
description of how to use the mass history models (analytic 
and semi-analytic) that were presented in Sections 2 and 4.4, 
respectively, see Appendix [Cl 
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Figure 8. Top panel: relation between formation redshift ( 2 — 2 ) 
and halo mass at 2 : = 0 (Mq). The different symbols correspond 
to median values, the error bars to la confidence limits and the 
grey area to the scatter. These were computed from the dark 
matter only simulations that assumed a WMAP3 (light blue line), 
WMAP5 (blue line), WMAPl (green line), WMAP9 (purple line) 
and Planck cosmology (dark blue line). The solid lines are not 
fits, but predictions of the Z —2 — Afhalo relation given by eq. (0. 
We also indicated the different values of the constant of propor¬ 
tionality A obtained by fitting the (p— 2 } — Pcrit(’2^— 2 ) relation. 
Bottom panel: halo mass history of a halo of 10^^ Mq at 2 : = 0 
from DMONLY simulations with different cosmologies. The grey 
curves show that as long as a suitable concentration-mass relation 
is assumed for the cosmology under consideration, eqs. ll37t . ll38ll 
and ilMli give a good estimate of the mass history curve. 


3.6.2 Mass definition dependence 

So far our calculations have been based on a halo mass de¬ 
fined as the mass of all matter within the radius r 2 oo at which 
the mean internal density (p)(< r 2 oo) is a factor of A = 200 
times the critical density of the universe, pcrit (from now 
on we denote this halo mass by M 200 ). In the literature a 


number of values ha ve been used for A. Some authors opt 
to use A = 200 fe.g.. I Jenkins et ^ 2001 ) or A = 200flm(2) 
(e.g. NFW), while others fe.g.. iBullock et al]l200l[ l choose 
A = Avir according to th e spherical virialization criterion of 
iBrvan fc Normaj] (Il998l l. These definitions can lead to size¬ 
able differences in c for a given halo and, as discussed, the 
differences are also cosmology-dependent. 

In this section we study how the structural properties 
and mass accretion histories depend on the adopted mass 
definition. We analyse halo mass histories of relaxed ha¬ 
los using three different halo mass definitions. First, we use 
M 200 . Secondly, we use Mmean, which is the mass within 
the radius rmean for which the mean internal density is 200 
times the mean background density. Finally, Mvir is the 
mass within the radius r^ir for which the mean internal 
density is Ayir times th e critical density as determined by 
iBrvan fc Normaj] lll998li . Note that halo masses and radii 
are determined using a sp herical overdensity ro utine within 
the SUBFIND algorithm (ISoringel ef^jlj JgOOl]) centred on 
the main subhalo of the FoF halos ( Davis et aklllOSSl i. We 
perform all calculations for the three different halo defini¬ 
tions, taking the halo centre to be the location of the parti¬ 
cle in the FoF group for which the gravitational potential is 
minimum. 

Eq. ([§1) shows that the formation redshift is directly 
proportional to the mean density within the scale radius 
((1 + Z- 2 )^ oc (p}(< ^■- 2 )), but the constant of proportion¬ 
ality depends on the mass definition that is adopted. There¬ 
fore, a change in the mass definition changes the formation 
time as 


(l + Zf)A, _ / (p)(< r-2)Ai 

(1 + z /) a 2 ~ ^-2)^2 


1/3 


(27) 


where A 1.2 refers to different overdensity criteria. That is, 
from eq. ([5]), (p)(< r- 2 ) changes according to the mass 
definition as (p)(< r_ 2 )A = A x pcrit.oCAU(l)/U(cA). Then 
(p)(< r_ 2 )Ai =200 refers to the mean internal density within 
r- 2 , obtained by defining the mean internal density at 
the virial radius to be 200pcrit,o. If we consider Ai = 200 
and A 2 = mean = 200S7m, and that there is a factor 
of 0.55 difference between the concentrations C200 and 
Cmean for a IO'^^Mq halo, then we obtain (using eq. I27II 
the relation (1 + Zf) 2 oo « (0.255n“^)^/^(l -I- 2:f)mean, 
where we have used the fact that 
(p)(< r_2)mean ~ (IIm/0.255)(p)(< r_2)20cEl. This im¬ 
plies that the change in the halo mass history due to 
different halo mass definitions is 


fogio 


M(z)200 

M (2^)mean 


Ioglo(l -I- 2)“200-«mean 
-l-logioe'^^™“'’”°“’^ (28) 

0.9560200 logio(l + 2 ) (29) 

-h[l - (0.255/fIn.)^/®] Iogio(e)/32oo2 
0.0543Z, (30) 


where in step (l2^ we replaced 0200 — Omean « 0.9560200, 


^ In this last step we used the approximation that (p)(< 2 )a = 

A X c^y(l)/y(c)pcrit,o — A X 0.643c^-^^Pcrit,o 
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Figure 9. Mass history of a IO^^Mq halo as a function of red- 
shift. The different coloured lines show the change in the mass 
history when different halo mass definitions are used. The green 
line shows the mass history of a halo of M 200 = IO^^Mq at 
2 = 0, whereas the dark blue (purple) line shows the mass history 
of a halo of Mmean = (Mvir = IO^^Mq) at 2 : = 0. The 

dashed lines show the mass history predicted by eqs. llWt . ll38t 
and (Enj. The difference lies in the formation redshift definition 
which is affected by the change in the mean inner density (see 
eq. [HJ. The value of (p)(< 'r_ 2 ) changes with the value of A we 
used in the definition of halo mass. The c — M relation correspond 
to the mass definition under consideration. 


which is valid for a halo, and a 200 = 0.2501 and 

P 200 = -0.8147. 

The difference in mass history given by eq. (O can 
be seen in Fig. [21 which shows how the halo mass his¬ 
tory is affected by the halo mass definition. The green 
line in Fig. [9] shows the mass history assuming the M 200 
mass definition. The purple line shows the Mvir defini¬ 
tion, and the dark blue line shows the Mmean definition. 
The different dashed lines correspond to the mass histories 
M{z) = M{z = 0)(1 -I- where the difference lies 

in the mass definition that changes the mean inner density 
and the concentration-mass relation (for a relaxed halo sam¬ 
ple). Duffy et al. (2008) studied how the halo mass definition 
changes the concentration-mass relation, and provided the 
parameters of the different c — M relations. They found that 
the concentration of a relaxed Mmean halo is 80% larger than 
the concentration of a relaxed M 200 halo. We adopt those 
fits in our calculations of the M{z) estimate and conclude 
that, as long as we use a concentration—mass relation that 
is consistent with the adopted halo mass definition, the ex¬ 
pressions (|S7|). (EHI and (ESI) accurately reproduce the halo 
mass history. 

The analytic estimate given by eq. (0 predicts that 
the difference in mass history due to the change in the mass 
definition (M 200 versus Mmean) is Alog;,o~ 0.05432. 
This can be seen in Fig. [S] where AlogjQM( 2 ) = 0.1086 
at 2 = 2 (Mmean (2 = 2) = and 

M2oo(2 = 2) = 10^-^i^Mq). 


4 COMPARISON BETWEEN SEMI-ANALYTIC 

AND ANALYTIC MODELS 

In this section we compare the semi-analytic model derived 
in Section rS.dl with the analytic model for halo mass history 
derived in Paper I. Note that while the semi-analytic model 
is obtained through fits to simulations, the analytic model 
was based on the extended Press-Schechter (EPS) theory 
without calibration against simulations. 

Fig. llOl shows a comparison between the models for var¬ 
ious halo masses (logiQ[Mo/MQ] = 5 — 15). As can be seen 
from the figure, the models mostly agree on the mass his¬ 
tories of halos with final masses between 10^'^ and IO^^Mq. 
However, there are a few factors of difference in the mass 
histories of larger and smaller halos, and the difference in¬ 
creases towards high redshift. We find that in the analytic 
model the halo mass decreases quite abruptly at high red¬ 
shift for halos with final masses > IO^^Mq. For instance, 
there is a factor of 9 difference at 2 = 5 between the models 
for a 10^® Mq halo. This difference is probably due to the 
progenitor definition. In the analytic model, the progenitor 
is defined as the halo with mass a factor of q lower (q ~ 4 
for Mo > IO^^Mq) at redshift 2 / ( 2 / ~ 0.9), whereas in the 
semi-analytic model, the progenitor is the halo that contains 
most of the 25 most bound dark matter particles in the fol¬ 
lowing snapshot. 

Fig.dSlalso shows that the semi-analytic model overpre¬ 
dicts the mass histories of halos with final masses < 10®Mq. 
This is expected because the parameters a and /3 in the semi- 
analytic model depend on the concentration-mass relation 
adopted. In this case we are using Duffy et al. (2008) rela¬ 
tion, which is calibrated in the mass range 10 ^'’ — IO^^Mq, 
for lower masses the con centration-mass rel ation deviates 
from a simple power law dLudlow et al.l[^14h . 

In Paper I we developed an analytic model of halo 
mass history, based on EPS theory, that only depends on 
the power spectrum of the primordial density perturbations. 
We found very good agreement between the halo mass his¬ 
tories predicted by o ur analytic model an d published fits 
to simulation results Jv ari den Boschl I 2 OO 2 I : iMcBride et al.l 
I 2 OO 9 I: I van den Bosch et al.ll2014h . In this work we have de¬ 
veloped a semi-analytic model that uses the functional form 
for the mass history motivated by EPS theory, and linked 
the mass history to halo structure through empirical rela¬ 
tions obtained from simulations. 

We now combine these two models (semi-analytic and 
analytic) to establish the physical link between a halo con¬ 
centration and the linear rms fluctuation of the primordial 
density field. From Fig. [10] we have found that there are a 
few factors of difference between the models. We now focus 
on the mass range 10^^ — 10^®Mq, where the factor of dif¬ 
ference is less than 1.5. We set the mass history curve to be 
the same in the two models, that is 


Af ( 2 ) Analytic — Af ( 2 ) Semi —analytic 


(31) 


for all redshifts. We then evaluate this equality at redshift 1 
and obtain. 


/(Mo)(o.92^U=o-0.3) = as(c)ln(2)-b/3s(c). (32) 

In this last equation, as(c) and /3s (c) are given by eqs. (1191) 
and m, respectively, and depend on concentration, D is 
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Figure 10. The top panel shows a comparison between the semi- 
analytic (solid lines) and the anal 3 d.ic model (dashed lines) for 
halo mass history. The bottom panel shows the residuals between 
the models. The different coloured lines correspond to the models 
for various halo masses in range log]^Q[Mo/MQ] = 5 — 15 in steps 
of 2. 


the linear growth factor and /(Mo) depends on the rms 
of the primordial density field, a. We approximate various 
terms in eq. ([3^ . including /(Mq) ^ 1.155(cr(Mo)^)^'^^^ and 
y(l)/y(c) 0.643c“°'^^, and obtain 

C = + 2.3, (33) 

which is suitable for a WMAP5 cosmology. Note that 
eq. (1331) is not a fit to any simulation data, it has been 
derived from eq. (1321) . Fig. 1111 shows the c — a relation at 
z = 0. In this figure we compare the predicted relation (solid 
line), as given by eq. (1331) (obtained by equaling the ana¬ 
lytic and semi-analytic models), with the simulation outputs 
(coloured symbols). The different symbols correspond to the 
median values of the relaxed sample of halos and the error 
bars to Icr confidence limits. The good agreement between 
the analytic prediction and the simulation outputs clearly 
shows that the halo mass accretion history is the physical 
connection in the c — cr correlation. 


5 SUMMARY AND CONCLUSION 

In this work we have demonstrated that there is an intrinsic 
relation between halo assembly history and inner halo struc¬ 
ture, and that the mass history is the physical connection 
between the inner halo structure and the power spectrum of 
initial density fluctuations. 

We examined the density profiles and mass growth 
histories of a large sample of halos and their progenitors 



-0.4 -0.2 0.0 0.2 0.4 

Log,oa“' 


Figure 11. Comparison between the a — c relation at 2 = 0 
predicted by the combination of the mass history models (solid 
line), and the simulation outputs (coloured symbols). 


within the OWLS simulations. We separated our halo sam¬ 
ple into a ‘relaxed’ sample, and a ‘complete’ sample that 
includes bot h relaxed and unrela xed halos. We confirmed 
the finding of iLudlow et al.l ll2013h that for relaxed halos the 
mean enclosed density within the NFW scale radius {r- 2 ), 
{p){< r- 2 ), is directly proportional to the critical density of 
the Universe at the formation redshift, z- 2 , dehned as the 
time at which the mass of the main progenitor equals the 
mass enclosed within the scale radius at z = 0, 

(p)(<r_ 2 ) = 900pcrit(2-2). (34) 

In the above relation, the value, 900±50, is obtained through 
hts to the simulation data (suitable for WMAP5 cosmology). 
We showed in Section [3.6.11 that expression (1341) provides a 
straightforward relation between formation time and con¬ 
centration, 

(1 -b z- 2 ) = [(2oo/9oo)c3y(i)/y(c) - ( 35 ) 

where y(c) = ln(l -|- c) — c/(l -|- c). The overall trend of 
decreasing formation redshift with decreasing concentration 
is evident for the main branches of the descendant halos. 
This implies that halos which assemble earlier have a higher 
concentration, because the density of the Universe at the 
formation time was larger. To relate the formation time to 
halo mass, we used the following concentration mass relation 
for relaxed halos 

c = 6.67(M/2 X (36) 

obtained bv iDuffv et al.l (l2008f) by fitting to the simulation 
data and suitable for the WMAP5 cosmology. We found 
that, on average, halo concentrations differ by a factor of 
1.16 between the relaxed and complete samples. The lower 
individual concentrations of unrelaxed halos (due to spuri¬ 
ous subhalos or ongoing mergers that do not result in an 
accurate fit for an NFW density profile) produce incorrect 
enclosed halo masses and therefore lower formation times 
(by a factor of 1.1). However, on average, the formation 



























The accretion history of dark halos II 15 


time-concentration relation does not change, thus indicat¬ 
ing that the halo mass history is not affected by the fact 
that a halo is out of equilibrium at a particular redshift. We 
used these findings to show that scatter in the halo mass 
history leads to scatter in the formation time {5z-2), and 
hence to scatter in the concentration—mass relation {5c). 

We found that formation time decreases with increasing 
mass (at a non-linear rate). This means that high-mass halos 
are still accreting mass rapidly in the present epoch, while 
low-mass halos typically accreted their mass early. Thus, the 
formation time—concentration relation provides the physical 
link between the halo mass history and internal structure. 
This result led us to provide a semi-analytic model for the 
halo mass history, which uses a direct, analytic correlation 
between the parameters a and /3 in the mass history leg. 1371) 
and concentration. 


M{z) = Mo{l + z)°‘e^\ (37) 

a = [ln(y(l)/y(c)) - , 92 - 2 ]/ln(l -I- 2 - 2 ), (38) 

13 = -3/(1+ 2 - 2 ), (39) 

where we obtained the constant value, —3.01 ± 0.08, in the 
last relation (eq. I39II by htting the halo mass history model 
(eq. 1371) to the simulation data. Then we obtained a semi- 
analytic model for the mass accretion history, that adopts 
the functional form, eq. (E3, and the parameters a and f3 
are now given by analytic relations that include numbers ob¬ 
tained from fits to simulation results. It is important to note 
that the semi-analytic model was derived assuming that the 
density profile of all halos is NFW. Interestingly, we found 
that the semi-analytic model describes with high-accuracy 
the mass histories of both relaxed and unrelaxed halo sam¬ 
ple, even though the NFW profile is not a good fit to the 
density profile of unrelaxed halos. 

We investigated how cosmology affects the semi- 
analytic model. We found that as long as a suitable 
concentration-mass relation and the value for the best-fitting 
parameter in the {p){< r- 2 ) — pcYit{z- 2 ) relation are as¬ 
sumed for the cosmology being considered, the semi-analytic 
models describe the mass histories with high accuracy. In ad¬ 
dition, we investigated how different mass definitions change 
the halo mass histories and we found that as long as we 
use a concentration-mass relation that is consistent with the 
adopted halo mass definition, the semi-analytic model accu¬ 
rately reproduces the halo mass history. 

In Paper I, we presented an analytic model for the halo 
mass history, based on extended Press-Schechter theory and 
not calibrated against simulations data. We compared the 
analytic model of Paper I with the semi-analytic model pre¬ 
sented here and found very good agreement in the mass 
range 10® — IO^^Mq. However, we found that the analytic 
model predicts larger masses at high redshift for halos with 
final masses > IO^^Mq, whereas the semi-analytic model 
overpredicts the mass history of low-mass halos (halos with 
final masses < 10®Mq). This is expected because the semi- 
analytic model depends on the adopted concentration-mass 
relation, which deviates from the assumed power law at low 
masses. The reader may find a step-by-step guide on how to 


implement the semi-analytic model in Appendix [C] as well 
as numerical routines onlin^^. 

Interestingly, by combining these two models (semi- 
analytic and analytic) we established the physical link be¬ 
tween a halo concentration and the initial density perturba¬ 
tion field, which explains the correlation between concentra¬ 
tion and rms fluctuation of the primordial density field, a 

(Fig.ini)- 

Finally, by differentiating eq. (EZl) we obtained the dark 
matter accretion rate, 

^ = 71.6MQyr-iMi2ho.7 

X [ — Q — /3(1 -|- 2)][f2m(l -|- 2 )® -|- (40) 

As the change in the a and jS parameters over halo masses is 
not significant when calculating accretion rates, we provided 
a mean accretion rate, obtained by averaging the parameters 
a and j3 over halo mass, 

(^> = 71.6MQyr-iMi2ho.7 

X[-0.24 + 0.75(1 + 2)][!:!m(l -f zf + IIa]^/^41) 

We then concluded that in order to predict halo mass 
growth, the concentration-mass relation from a relaxed sam¬ 
ple should be used. 

Putting the pieces together, we addressed the question 
of how the structure of halos depends on the primordial 
density perturbation field. We found that concentration is 
the link between the halo mass profile and the halo mass 
history (and that one can be determined from the other). 
We also found that the ‘shape’ of the halo mass history is 
given by the linear growth factor and linear power spectrum 
of density fluctuations. Therefore, we concluded that halo 
concentrations are directly connected to the initial density 
perturbation field. 

In a forthcoming paper (Paper III) we will combine the 
analytic model presented in Paper I and semi-analytic model 
presented here to predict the concentration-mass relation. 
We will investigate its evolution. We will show that extrap¬ 
olations to low masses of power-law fits to simulation results 
are highly inadequate, and will investigate whether linear 
(p)(< r- 2 ) — Pcrit{z- 2 ) relation holds at redshifts other than 
0 . 
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tree. The different mass histories are calculated by following 
the merger trees of a given sample of halos. At each redshift 
the mass histories are computed by calculating the median 
mass value, determined by bootstrap resampling the halos, 
from the merger tree. Along with the median value, the la 
confidence interval is recorded. 


APPENDIX A: SIMULATIONS 


The DMONLY simulations contain only dark matter parti¬ 
cles, which provide us with a useful baseline model for testing 
the impact of baryonic physics on halo mass histories and 
mass profiles, when comparing it with the REF simulation. 
The REF simulation includes sub-gr id recipes for star forma¬ 
tion (ISchave fc Dalla Vecchia[ j2008[l. radiat ive (metal-line) 


cooling and heating i Wiersma et al.l 12009811 . stellar evolu- 


tion, mass loss from massive stars and chemical enrichment 


(IWiersma et al. 

l2009bri. and a kinetic prescription for super- 

nova feedback ( 

Dalla Vecchia & Schavell200^. In this work 


we use hydrodynamical simulations (along with DMONLY) 
to test the effect of baryons on halo mass histories and com¬ 
pute the gas accretion rate. We do not test other fe e dback 
schemes because it was shown by Ivan de Voort et al.l (I 2 OIIII 
that the halo mass accretion is robust to variations in feed¬ 
back. 


A1 Construction of halo merger trees 

The first step towards studying the mass assembly history 
of halos is to identify gravitationally bound structures and 
build halo merger trees. We begin by selec ting the largest 
halo in each FoF group (I Davis et al.l Il985ll (i.e. the main 
subhalo of FoF groups that is not embedded inside larger 
halos). Halo virial masses and radii are determined using 
a sphe rical overdensity routine within the S UBFIND algo¬ 
rithm (ISpringel et al.l[200ll : lDolag et al.ll2009l i centred on the 
main subhalo of FoF halos. Therefore, we define halo masses 
as all matter within the radius r 2 oo for which the mean in¬ 
ternal density is 200 times the critical density. Throughout 
this work, we study the accretion history of the largest ha¬ 
los in each FoF group. Subhalos, defined as bound structures 
that reside within the virial radius of the largest ‘host’ halo, 
show distinct mass histories. The structures of subhalos are 
strongly affected by the potential of their host halos, as seen 
for example in the cessation of mass ac cretion on to subh a¬ 
los residing in dense env ironments fsee lWang et al.ll2009l or 
iLacerna fc Padillall201lh . Consequently, the masses of sub¬ 
halos do not follow the mass history of their host halos. 

The merger trees of the largest halos are then built as 
follows. First, at each output redshift (snapshot), we select 
‘resolved’ halos that contain more than 300 dark matter par¬ 
ticles. We refer to these resolved halos as ‘descendants’. We 
then link each descendant with a unique ‘progenitor’ at the 
previous output redshift. This is nontrivial due to halo frag¬ 
mentation: subhalos of a progenitor halo may have descen¬ 
dants that reside in more than one halo. The fragmentation 
can be spurious or due to a physical unbinding event. To 
correct this, we link the descendant to the progenitor that 
contains the majority of the descendant’s 25 most bound 
dark matter particles. Therefore, the main progenitor of a 
given dark matter halo is found by tracing backwards in time 
to the most massive halo along the main branch of its merger 


APPENDIX B: ANALYSIS OF SCATTER 

BO. 1 Scatter in formation times and concentrations 

We now analyse the scatter in the formation time—mass 
relation and show it relates to the scatter in the 
concentration—mass relation. So far, we have shown that 
the formation time is related to halo concentration through 
the mean inner density ((p)(< r_ 2 ))—critical density 

(pcrit( 2 - 2 )) relation plotted in Fig. [T] Through first order 
error propagation, we look for the corresponding scatter in 
the formation time. 
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^ J [ D ,„(1 + z_2)® + Ba] 



3Dm(l + Z— 2 ) 2 — 2 


5z.2\ 

(p)(< r_ 2 ) /Jc\ 

U-2 j 

900po V c / 


(Bl) 

(B2) 

(B3) 
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(1 -I- cY Y{c) 

where we used S{c^/Y{c)) = ZYSc/Y{ c) — c^Sc/\Y{cY{ 1 + 
cY]. Eq. (IB4II relates the scatter in formation time 
{\dz- 2 \lZ-D to the scatter in the concentration (|Jc|/c). 

The grey shaded areas in the panels in Fig. [3] show the 
scatter in 2 _ 2 , while the panels in Fig. lAll show the scatter 
in the c — M relation (left panel) and in the Z -2 — M relation 
(right panel). The grey contours in Fig. I All enclose 68% of 
the distribution while the individual points show the remain¬ 
ing 32%. The black (green) solid line shows the mean scatter 
in the formation time per halo mass bin for the relaxed (com¬ 
plete) sample. The presence of unrelaxed halos does not have 
any significant effect on either the scatter in the formation 
time or the mass histories. The average scatter in formation 
time is t\&z- 2 \lz- 2 ) = 0.324 {{\5z-2\/z- 2 ) — 0.356) for the 
relaxed (complete) sample. 

The left panel of Fig. lAll shows that the average scatter 
in the concentration—mass relation is (|Jc|/c) = 0.257 for 
the full sample and (|Jc|/c) = 0.218 for t he relaxed sample . 
In agreement with previous work (see e.g. lNeto et al.ll2007h . 
the scatter in the concentration of the relaxed halo sample is 
lower than the scatter of the full sample. The extra scatter 
in the full sample is produced by the deviation of the density 
profiles from the NFW form for halos experiencing ongoing 
mergers and for artificially linked halos. 

Assuming {\Sz- 2 \/Z- 2 ) ~ 0.324, {\5c\/c) can be ob¬ 
tained as a function of halo mass by applying eq. (fBil) . This 
analytic estimate is plotted in the left panel for the relaxed 
sample (red dashed line). We find very good agreement be¬ 
tween the scatter in concentration from eq. (IB4II . and the 
median value plotted in black for the relaxed sample, and in 
green for the complete sample. Therefore, we conclude that 
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Figure Al. Scatter in the concentration—mass relation (left panel) and the formation redshift—mass relation (right panel). Left panel: 
Y-axis shows the difference between the concentration predicted by Duffy et al. (2008) (see eq. [3|) and the actual concentration from 
the simulation output. This difference is divided by the concentration from the simulation output and plotted against virial halo mass. 
The different colours of the points indicate that the concentration outputs were obtained from DMONLY_WMAP5 simulations with 
different box sizes. The grey contours enclose 68% of the distribution while the individual points show the remaining 32%. The black 
(green) solid line shows the mean relative scatter in the concentration-mass relation per halo mass bin of the relaxed (complete) sample. 
The red dashed line is an analytic estimate of the scatter obtained by propagating the scatter in the formation redshift—mass relation 
to the concentration—mass relation (eq. iBit . Right panel: same as the left panel but the scatter is obtained from the difference in the 
formation redshift predicted from eq. m and the simulation output. The black (green) solid line shows the median value in the scatter 
as a function of mass of the relaxed (complete) sample. 


the scatter in formation time determines the scatter in the 
concentration. However, at higher mass es and redshifts the 
fraction of relaxed halos decreases (e.g. iLudlow et aLllioi^ ) 
and it has been fonnd that the concentration-mass relation 
of the c omplete halo sample exhibits a strong flattening and 
upturn llKlvpin et al.l[20ll] : IPrada et al.ll2012) . As a result, 
at high masses and redshifts the scatter in the concentration 
will probably depend on other variables besides the scatter 
in formation time. 

In the following subsection we find that the scatter in 
the accretion history determines the scatter in the formation 
time. 


BO. 2 Scatter in halo mass histories 

In this section we analyse the scatter found when computing 
mass histories from the simulation outputs. We analytically 
estimate the scatter in the mass history due to both the scat¬ 
ter in the concentration and formation times (estimated in 
Section rB0.il) . We then compare this to the scatter obtained 
from the simulation outputs. 

To compute the scatter in the mass history we perform 
a first order error propagation in M{z) (eg. HOI) . 

SM{z)/M{z) = 5a\n{l + z) + zSP, (B5) 

where 5a {SP) is the scatter in a (P) due to the scatter in 
c and Z- 2 . From eq. m we first compute Sa due to the 
scatter in Z-2, 


= S{-Pz-2l ln(l -I- Z-2)) 


— Z- 2 SP — PSz -2 
ln(l -I- Z- 2 ) 


Pz- 25 fn{l + Z-2)] ^ 


Z-2PA + z-2\n{l + Z-2)) ^ 5 z- 2 \ 
(1 + Z-2) la(l + Z-2) \ Z-2 )’ 


(B6) 


where in the last line we used SP = ,, 2. ) which 

(l-fZ-2) 2-2 J ’ 

follows from eq. m- Similarly, we calculate the scatter in 
a due to the scatter in c, 


5a = 5(ln(y(l)/y(c))/ln(H-z_2) 

5Y{c) 1 

y(c) ln(l -I- Z-2) 


(1 -I- c)2y(c) ln(l + Z-2) \ c J ’ 


(B7) 


where in the last step we used SY{c) = (^)- Fi“ 

nally, we substitute Sp = ) and eqs. (IB 6 II and 

(TBTI) into eq. HSl). From Section IBP. II we adopt the av¬ 
erage scatter in Z -2 and c for the complete sample, that 
is (|5c|/c) = 0.25 and (5z_2/z-2) = 0.35, and compute 
{\dM{z)/M{z)\) as a function of redshift. Fig. IBII shows 
the scatter measured from the simulation outputs (coloured 
lines) and from the analytic estimates (grey dashed lines). 
The different coloured lines correspond to the median val¬ 
ues of the scatter from simulations with different box sizes. 
We calculate these by averaging the difference between the 
mass history predicted by eqs. (fTO]). ([Till and (I 2 III (for a given 
halo with mass M 200 at 2 ; = 0), and each M{z) given by the 
merger trees from the complete halo sample. We find very 
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Figure Bl. Mean scatter in the halo mass history against redshift 
for different halo masses M2oo- The y-axis shows the mean value 
of the difference between the mass history (M^t) predicted by 
eqs. (llOll . (I14II and (I21II . and the mass history (Msim) from the 
simulation output. The difference is divided by the Msim from 
the simulation output. The different coloured lines correspond to 
the different DMONLY_WMAP5 simulations as indicated in the 
legends. The grey lines show the analytic estimates of the mass 
history curves given by eqs. llBSll . llB6ll and 


good agreement indicating that the scatter in the concen¬ 
tration comes from the scatter in the formation time, which 
in turn comes from the scatter in the halo mass history. 


APPENDIX C: STEP-BY-STEP GUIDE TO 
COMPUTE HALO MASS HISTORIES 

CO. 3 Semi-analytical model 

This Appendix provides a step-by-step procedure that out¬ 
lines how to calculate the halo mass histories using the nu¬ 
merical model presented in Section 4.4: 

(i) First assume a cosmology and choose a concentration- 
mass relation from the literature. For instance, the Duffy 
et al. (2008) relation, c = 6.67(Mo/2 x 

is suit able for the WMAP5 cosmology, whereas iNeto et al.l 
ll2007tl is suitable for WMAPl: 

(ii) Calculate the formation time, 

( 200 c{MofY{l) 

" ^AcosmoDmY(c(Mo)) ^ ’ 

Note that the value of Acosmo in the above equation is 
cosmology dependent. Awmaps = 900 is suitable for the 
WMAP5 cosmology. In this work we obtained Awmapi = 
787, Awmaps = 850, Awmap9 = 820 and Apianck = 798. 

(hi) Calculate the parameters a and /3, 

Q = [ln(y(l)/y(c))-/3z_2]/ln(l + 2 _ 2 ), (C2) 

p = -3/(1-b^-2). (C3) 

(iv) Finally, the mass history can be calculated as follows, 


M(2) = Mo(l-bz)“e^y (C4) 

The above model is suitable for any cosmology (as long 
as the concentration-mass relation and the value of Acosmo 
correspond to the desire cosmology) and is valid over the 
halo mass range for which the concentration-mass and the 
Z -2 — Mo relations, obtained from simulations, are valid (e.g. 
10^° — IO^^Mq for Duffy et al. 2008). 














